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The main contribution of this paper is to derive an explicit expression for the fundamental pre¬ 
cision bound, the Holevo bound, for estimating any two-parameter family of qubit mixed-states 
in terms of quantum versions of Fisher information. The obtained formula depends solely on the 
symmetric logarithmic derivative (SLD), the right logarithmic derivative (RLD) Fisher information, 
and a given weight matrix. This result immediately provides necessary and sufficient conditions 
for the following two important classes of quantum statistical models; the Holevo bound coincides 
with the SLD Cramer-Rao bound and it does with the RLD Cramer-Rao bound. One of the im¬ 
portant results of this paper is that a general model other than these two special cases exhibits an 
unexpected property: The structure of the Holevo bound changes smoothly when the weight matrix 
varies. In particular, it always coincides with the RLD Cramer-Rao bound for a certain choice of 
the weight matrix. Several examples illustrate these findings. 


I. INTRODUCTION 

One of the fundamental questions in quantum statisti¬ 
cal inference problem is to establish the ultimate precision 
bound for a given quantum statistical model allowed by 
the laws of statistics and quantum theory. Mainly due 
to the non-commutativity of operators and nontrivial op¬ 
timization over all possible measurements, this question 
still remains open in full generality. This is very much in 
contrast to the classical case where the precision bounds 
are obtained in terms of information quantities for vari¬ 
ous statistical inference problems. 

The problem of point estimation of quantum paramet¬ 
ric models is of fundamental importance among various 
quantum statistical inference problems. This problem 
was initiated by Helstrom in the 1960s and he devised a 
method to translate the well-known strategies developed 
in classical statistics into the quantum case [1]. A quan¬ 
tum version of Fisher information was successfully intro¬ 
duced and the corresponding precision bound, a quan¬ 
tum version of Cramer-Rao (CR) bound, was derived. 
It turned out, however, that the obtained bound is not 
generally achievable except for trivial cases. 

A clear distinction regarding the quantum parameter 
estimation problem arises when exploring possible esti¬ 
mation strategies since there is no measurement degrees 
of freedom in the classical estimation problem. Consider 
n identical copies of a given quantum state and we are 
allowed to perform any kinds of quantum measurements 
according to quantum theory. A natural question is then 
to ask how much can one improve estimation errors by 
measurements jointly performed on the n copies when 
compared to the case by those individually performed on 
each quantum state. The former class of measurements 
is called collective or joint and the latter is referred to as 
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separable in literature. It is clear that the class of collec¬ 
tive measurements includes separable ones and one ex¬ 
pects that collective measurements should be more pow¬ 
erful than separable ones in general. Since one cannot 
do better than the best collective measurement, the ul¬ 
timate precision bound is the one that is asymptotically 
achieved by a sequence of the best collective measure¬ 
ments as the number of copies tends to infinity. This 
fundamental question has been addressed by several au¬ 
thors before [2-14]. 

It was Holevo who developed parameter estimation the¬ 
ory of quantum states by departing from a direct analogy 
to classical statistics. He proposed a bound, known as the 
Holevo bound, in the 1970s aiming to derive the funda¬ 
mental precision limit for quantum parameter estimation 
problem, see Ch. 6 of his book [15]. At that time, it 
was not entirely clear whether or not this bound is a re¬ 
ally tight one, i.e., the asymptotic achievability by some 
sequence of measurements. Over the last decade, there 
have been several important progress on asymptotic anal¬ 
ysis of quantum parameter estimation theory revealing 
that the Holevo bound is indeed the best asymptotically 
achievable bound under certain conditions [10-13]. These 
results confirm that the Holevo bound plays a pivotal role 
in the asymptotic theory of quantum parameter estima¬ 
tion problem. Despite the fact that we now have the 
fundamental precision bound, the Holevo bound has a 
major drawback: It is not an explicit form in terms of a 
given model, but rather it is written as an optimization of 
a certain nontrivial function. Therefore, unlike the clas¬ 
sical case, where the Fisher information can be directly 
calculated from a given statistical model, the structure 
of this bound is not transparent in terms of the model 
under consideration. 

Having said the above introductory remarks, we wish 
to gain a deeper insight into the structure of the Holevo 
bound reflecting statistical properties of a given model. 
To make progress along this line of thoughts, we take 
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the simplest quantum parametric model, a general qubit 
model, and analyze its Holevo bound in detail. Since ex¬ 
plicit formulas for the Holevo bound for mixed-state mod¬ 
els with one and three parameters and pure-state models 
are known in literature, the case of two-parameter qubit 
model is the only one left to be solved. The main contri¬ 
bution of this paper is to derive an explicit expression for 
the Holevo bound for any two-parameter qubit model of 
mixed-states without referring to a specific parametriza- 
tion of the model. Remarkably, the obtained formula 
depends solely on a given weight matrix and three previ¬ 
ously known bounds: The symmetric logarithmic deriva¬ 
tive (SLD) CR bound, the right logarithmic derivative 
(RLD) CR bound, and the bound for D-invariant models. 
This result immediately provides necessary and sufficient 
conditions for the two important cases. One is when the 
Holevo bound coincides with the RLD CR bound and 
the other is when it does with the SLD CR bound. We 
also show that a general model other than these two spe¬ 
cial cases exhibits an unexpected property, that is, the 
structure of the Holevo bound changes smoothly when 
the weight matrix varies. We note that similar transi¬ 
tion has been obtained by others [10, 13] for a specific 
parametrization of two-parameter qubit model. Here we 
emphasize that our result is most general and is expressed 
in terms only of the weight matrix and two quantum 
Fisher information. 

The main result of this paper is summarized in the 
following theorem (The detail of these quantities will be 
given later): Consider a two-parameter qubit model of 
mixed states, which changes smoothly about variation of 
the parameter. Denote the SLD and RLD Fisher infor¬ 
mation matrix by Gg and Gg , respectively, and define the 
SLD and RLD CR bounds by 

C$[W]=Tt{WG^}, (1) 

C?[W\ = Tr | IT Re Gg 1 } + TrAbs { WIm Gg 1 } , (2) 


qubit model under the regularity conditions is 


C?[W} = 


C$[W] if G e fl [W] > 


CZ[W] + Ci[W\ 


G,f [W] + S'efW] otherwise 


( 5 ) 


where the function SgfW], which is nonnegative, is de¬ 
fined by 


Sg[W} := 


[■\(Cg[W}+C§[W])-Cf[W ]] 2 
Of [W] - C«[W] 


( 6 ) 


Note that the condition Cf < (Gg +Cg)/2 implies Cg — 
C*> 0 so that Sg [W] is well defined. (See the discussion 
after Eq. (75).) 

The above main result 1.1 sheds several new insights 
on the quantum parameter estimation problems. First 
note that the form of the Holevo bound changes accord¬ 
ing to the choice of weight matrices. This kind of transi¬ 
tion phenomenon has never occurred in the classical case. 
Second surprise is to observe the appearance of the RLD 
CR bound in the generic two-parameter estimation prob¬ 
lem. As we will provide in the next section, the RLD 
CR bound has been shown to be important for a special 
class of statistical models, known as a D-invariant model. 
Here, we explicitly show that it also plays a major role for 
non D-invariant models. Last, in many of previous stud¬ 
ies parameter estimation problems, the precision bound 
is either expressed in terms of the SLD or RLD Fisher 
information, but the second case of the above expression 
(5) depends both on the SLD and RLD Fisher informa¬ 
tion. To see this explicitly, we can rewrite it as 


C^[W]+S e [W]=Tr{WGg 1 } 

]_ (TrAbs { WIm G^ 1 }) 

4 TrjwtG^-ReG,- 1 )}' 


( 7 ) 


respectively. Here, W is a given 2x2 positive definite 
matrix and is called a weight matrix (TrAbs {} is defined 
after Eq. (22).). Introduce another quantity by 

Gf [W\ := Tr {WRe Zg) + TrAbs {WIm Zg} , (3) 

where Zg is a 2 x 2 hermite matrix defined by 

z e : = [ z e]i, je {i,2y z e ;= tr (po L e L e) ■ ( 4 ) 

Here Lg are a linear combination of the SLD operators; 
Lg = ^2j 9g Lg,j with g l g J denoting the (i,j) component 
of the inverse of SLD Fisher information matrix and Lg t i 
SLD operators. With these definitions, we obtain the 
following result: 

Theorem 1.1 The Holevo bound for any two-parameter 


All these findings will be discussed in details together 
with examples. 

The rest of this paper continues as follows. Section II 
provides definitions and some of known results for pa¬ 
rameter estimation theory within the asymptotically un¬ 
biased setting. In Sec. Ill, a useful tool based on the 
Bloch vector is introduced and then the above main the¬ 
orem is proved. Discussions on the main theorem are 
presented in Sec. IV. Section V gives several examples to 
illustrate findings of this paper. Concluding remarks are 
listed in Sec. VI. Most of the proofs for lemmas are de¬ 
ferred to Appendix A. Supplemental materials are given 
in Appendix B. 

II. PRELIMINARIES 

In this section, we establish definitions and notations 
used in this paper. We then list several known results 
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regarding the Holevo bound to make the paper self- 
contained. 


A. Definitions 

Consider a d-dimensional Hilbert space 'H (d < oo) and 
a fc-parameter family of quantum states pg on it: 

M :={p e \9={9 1 ,9 2 ,...,9 k ) G 0cK fc }, (8) 

where 0 is an open subset of fc-dimensional Euclidean 
space. The family of states M. is called a quantum sta¬ 
tistical model or simply a model. The model discussed 
throughout the paper is assumed to satisfy certain regu¬ 
larity conditions for the mathematical reasons [16]. For 
our purpose, the relevant regularity conditions are: i) 
The state pg is faithful, i.e., pg is strictly positive, ii) It 
is differentiable with respect to these parameters 6 suf¬ 
ficiently many times, iii) The partial derivatives of the 
state dpg/dd 1 are all linearly independent. In the rest of 
this paper, the regularity conditions above are taken for 
granted unless otherwise stated. 

For a given quantum state p G S(H), we define the 
SLD and RLD inner products by 

(X, Y) p := itr (p(YX* + X*Y)) = Re {tr {pYX*)}, 

(X,Y)+ :=tv(pYX*), (9) 

respectively, for any (bounded) linear operators X, Y on 

'H. Here, * denotes the hermite conjugation. Given a 
^-parameter model A4 = {/90|0G0CR fc }, the SLD op¬ 
erators Lg i and RLD operators Lg t are formally defined 
by the solutions to the operator equations: 

d 1 

-qqP® = -j\Pe L e,i + Lg,ipe), 

d 

— po = pgLg^i. (10) 

The SLD Fisher information matrix is defined by [1] 

Gg := (-*-■*■) 

90,ij • (Lg,i, Lg j)p Q — tl' ^pg—(Lg,Lgj T LgjLg^j^ , 
and the RLD Fisher information is [15, 17] 

Gg := \i9o,ij]i t j e {i t 2,...,k} 

9o,ij := {Lg,i, Lg j)~^ e = tr ^pgLg jLg^j . (12) 

Define the following linear combinations of the SLD 
and RLD operators: 

k k 

Lg ■■= and Lg := ^(Gg^Lgj. 

1=1 l=i 


From these definitions, the following orthogonality con¬ 
ditions hold. 

(Li, Lg tj ) p , = 6), and ( U e , L ej )+ = 5). (13) 

These operators with upper indices are referred to as the 
SLD and RLD dual operators, respectively. 

Consider nth i.i.d. extension of a given state and we 
define nth extended model by 

M (n) :={pf n \9e 0cK fc }. (14) 

The main objective of quantum statisticians is to perform 
a measurement on the n tensor state pf n and then to 
make an estimate for the value of the parameter 9 based 
on the measurement outcomes. Here measurements are 
described mathematically by a positive operator-valued 
measure (POVM) and is denoted as n( n '. An estimator, 
which is a purely classical data processing, is a (measur¬ 
able) function taking values on 0 and is denoted as 9 n . 
They are 

= {n<>) | x g x n , > o, ^ njw = /® n }, 

xex n 

= (ole 2 n ,...,o k n ): x n ^e, 

where L is the identity operator on T-L and we assume that 
POVMs consist of discrete measurement outcomes. For 
continuous POVMs, we replace the summation by an in¬ 
tegration. A pair (n^, 9 n ) is called a quantum estimator 
or simply an estimator when is clear from the context and 
is denoted by ll^ [18]. 

The performance of a particular estimator can be com¬ 
pared to others based on a given figure of merit and then 
one can seek the “best” estimator accordingly. As there is 
no universally accepted figure of merit, one should care¬ 
fully adopt a reasonable one depending upon a given sit¬ 
uation. For example, a specific prior distribution for the 
parameter 9 is known, the Bayesian criterion would be 
meaningful to find the best Bayesian estimator. If one 
wishes to avoid bad performance of estimators, the min- 
max criterion provides an optimal one that suppresses 
such cases. In this paper, we are interested in analyz¬ 
ing estimation errors at specific point 9 G 0, that is, the 
pointwise estimation setting. For a given model (14) and 
an estimator IT”) = (n*”\ 9 n ), we define a bias at a point 
9 as 

«4 n) [n (n) ] := £ (9 n (x) - 9)p ( e n \x) = E ^[9 n ] - 9, 

xeXn 

with p^\x) := tr , (15) 

where E^fA*^] denotes the expectation value of a ran¬ 
dom variable X^ with respect to the probability dis¬ 
tribution pg l \ Note that the bias [£[(”)] is a k - 
dimensional real vector. An important class of estimators 
when estimating the specific point of the model is the lo¬ 
cally unbiased estimatior. This is to restrict estimators 
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such that the bias vanishes at the true point 9 up to the 
first order in the Taylor expansion. Mathematically, an 
estimator IT n ) is called locally unbiased at 9 if 

b\ [n (n) ] = 0 and A E («) { 0 O n] = S {, (16) 

hold at 9 £ 0. It is known that the Quantum CR bounds 
hold for any locally unbiased estimator [1, 15]. 

Upon analyzing performance of estimators within the 
asymptotic regime, we should impose some conditions 
that restrict the class of estimators. In statistics, a se¬ 
quence of an estimator is said (weakly) consistent, if it 
converges to the true value in probability for every value 
9 £ 0, i.e., 

Ve > 0 lim Pr{|0„ — 9\ > e} = 0, 


holds for all 9 £ 0. In this expression, |u| = (Xw?) 1 ^ 2 
denotes the standard Euclidean norm and the right hand 
side means that error probability can be made arbitrary 
small. As a good estimator must converge to the true 
value as n goes to infinity, it is reasonable to look for the 
class of consistent estimators in quantum parameter es¬ 
timation as well. In classical statistics, this condition of 
consistency alone turns out to be weak in order to exclude 
artificial estimators. There are several approaches to han¬ 
dle these problems in the classical case [19]. Rather than 
going into mathematical discussions, we simply look for 
the following class of estimators to avoid such a situation. 
A sequence of estimators is called asymptoti¬ 

cally unbiased if it satisfies 

lim yn^ n) [f[ (n) ] = 0, I™ [0 3 n ] = 8 ^, (17) 

for all i,j £ {1, 2,..., fc} and for all 9 £ 0. That is to 
require the locally unbiased condition (16) in the n —> oo 
limit. 

To quantify estimation errors of a given estimator, we 
consider the mean-square error (MSE) matrix defined by 


yin) 


V 


ij 

6,n 


- [<(»"’]] i j6(I 2 . 

[n (n) ] := ]T (%(x) - 9 i )(§i(x) - t 5»>< n) 


(x). 


x&X n 


By definition, the MSE matrix is a k x k real symmetric 
matrix and it is straightforward to show that it is non¬ 
negative. As stated in the introduction, we wish to find 
the best precision bound allowed by the laws of quantum 
theory and statistics, which is achievable in the n —> oo 
limit. In the classical case, one can directly minimize 
the MSE matrix as a matrix inequality over the class of 
asymptotically unbiased estimators and to find the lowest 
MSE error achievable as n —> oo. This line of approach 
does not work in the quantum case. One way to tackle 
this question is to deal with a weighted trace of the MSE 
matrix, which is a scalar quantity, and it is defined by 

MSE^ffl^lW] := TV |iuu e (n) [n (n) ]|. 


Here the matrix W is called a weight matrix and can be 
chosen arbitrary as long as it is real and strictly positive. 
Since the weight matrix is one of the important ingredient 
for our discussion, let us denote the set of all possible 
weight matrices by 

W := {W £ R kxk | W > 0}. (18) 

By changing the weight matrix, one can explore trade-off 
relations in estimating different parameter components 
9 l . We note that the weight matrix can depend on the 
value of the estimation parameter 9 as well. For example, 
it can be chosen as the SLD Fisher information matrix. 

Defining these terminologies, we now state the prob¬ 
lem: For a specific point of a given i.i.d. model Xi^ n \ 
what is the best sequence of estimators {IT”)}]^ and 
the minimum value of the weighted trace of MSE? To 
put it differently, one wishes to find the optimal sequence 
of estimators that minimum of the first order coefficient 
Ce[IU] in the large n expansion: 

MSE^ n) [n (n) |W] ~ 2-1 + 0{n~ 2 ), (19) 

i.e., the fastest decaying rate for the MSE. Mathemati¬ 
cally, we define the CR type bound for the MSE by the 
following optimization problem: 

Cg[W\:= inf { lim sup n MSE^ } [ft^ \W ]}, 
{n<»}:a.u. n-> oo 

where the infimum is taken over all possible sequences 
of estimators that is asymptotically unbiased 

(a.u.). Note that this bound depends both on the weight 
matrix W and the model pg at 9. The symbol 9 appearing 
in the bound C^flU] represents the model pg at 9. Unlike 
the Bayesian or the min-max settings mentioned before, 
we are interested in understanding statistical properties 
of a given parametric model. This would be important 
in particular study of quantum states from geometrical 
point of view [20]. 


B. The Holevo bound 

To define the Holevo bound, we need some defini¬ 
tions first. For a given quantum statistical model on 
H, denote a k array of hermite operators on H by 
X = (. X\X 2 ,...,X k ), X 1 £ C H (U), he, (X 1 )* = X\ 
for alii = 1,2,..., k, and define the set Xg by 

:= {X | \/i X 1 £ £ h (U),Vi tr (p e X l ) = 0, 

= ^ } - (20) 

The Holevo function [21] in the quantum estimation the¬ 
ory is defined by 

hg\X\W] := Tr jlURe Zg[X] j + TrAbs j Wlm Zg[X] j , 

(21) 
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where the k x k hermite matrix Zg[X] is 

Z e [X}:=[tr(p e X*X% Je{1 _ k} , (22) 

and TrAbsA denotes the sum of the absolute values of A j 
with X = Tdiag(Ai, A 2 ,..., A lrl )T^ 1 for some invertible 
matrix T. We note the following relation also holds for 
any anti-symmetric operator X: 

TrAbs {WX} = |A*| = Tr ||W 1/2 AW 1/2 ||, (23) 

i 

where |X| = \[X*X denotes the absolute value of a linear 
operator X. 

The Holevo bound is defined through the following op¬ 
timization: 

Cg[W] := min hg[X\W}. (24) 

XGX e 

The derivation of the above optimization is well summa¬ 
rized in Hayaslri and Matsumoto [10]. Holevo showed 
that this quantity is a bound for the MSE for estimating 
a single copy of the given state under the locally unbiased 
condition [15]: 

MSE^plW] > Cf [W], (25) 

holds for any locally unbiased estimator n. The nontriv¬ 
ial property of the Holevo bound is the additivity [10]: 

C?[W,pf n ]=n- 1 C?[W,p e ], (26) 

where the notation Cg [W, pf n ] represents the Holevo 
bound about nth extended model. 

The following theorem establishes that the Holevo 
bound is the solution to the problem of our interest: 

Theorem II. 1 For a given model satisfying the regular¬ 
ity conditions, Ce[W] = C^\W\ holds for all weight ma¬ 
trices. 

There exist several different approaches upon proving the 
above theorem. Hayashi and Matsumoto [10] proved the 
case for a full qubit model first. Gu^a and Kahn [11] 
introduced a different tool based on (strong) quantum 
local asymptotic normality to prove the qubit case. This 
was further generalized to full models on any finite di¬ 
mensional Hilbert space [12]. However, all these proofs 
depend on a specific parametrization of quantum states. 
More general proof has been recently established by Ya- 
magata, Fujiwara, and Gill [13]. 

This theorem implies that if we choose an optimal se¬ 
quence of estimators, the MSE behaves as 

Tr {H% (n) [fl (n) ]} ~ Cf> ^ + 0(n“ 2 ), (27) 

for sufficiently large n. That is the Holevo bound is the 
fastest decaying rate for the MSE. 


Although the Holevo bound stands as an important 
cornerstone to set the fundamental precision bound, the 
definition (24) contains a nontrivial optimization. The 
main motivation of our work, as stated in the introduc¬ 
tion, is to perform this optimization explicitly for any 
given model for qubit case. The result shows several 
nontrivial aspects of parameter estimation in quantum 
domain. Before going to present our result, we summa¬ 
rize several known results. 


C. Holevo bound for one-parameter and 
D-invariant models 

In this subsection, we consider two special cases where 
analytical forms of the Holevo bound are known. 

For a given k- parameter model on the Hilbert space Li, 
let us denote SLD and RLD Fisher information matrices 
by Gg and Gg, respectively, Eqs. (11, 12). Define the 
SLD and RLD CR bounds by 

C^m^TriWGg 1 }, (28) 

Cf[W] := Tr {wReGg 1 } + TrAbs {wimGff 1 } , (29) 

respectively. Throughout the paper, we use the notation 
ReG^ 1 = RejGg 1 } (IinG^ 1 = ImlG^ 1 }) representing 
the real (imaginary) part of the inverse matrix of the RLD 
Fisher information matrix. The well-known fact is that 
the SLD and RLD CR bounds cannot be better than the 
Holevo bound: 

Lemma II. 2 For a given model satisfying the regularity 
conditions, the Holevo bound is more informative than 
the SLD and the RLD CR bound, i.e., Cf\W\ > Cg\W\ 
and Cf[W] > G^[LF] hold for an arbitrary weight matrix 
W. 

Proof can be found in the original work by Holevo that 
is summarized in his book [15]. More compact proof was 
stated by Nagaoka [21]. See also Hayashi and Matsumoto 
[ 10 ]. 

1. One-parameter model 

When the number of parameters is one, the problem 
can be reduced significantly. In this case, there cannot 
be any imaginary part for the matrix (22) and thus the 
minimization is reduced to minimizing the MSE itself. 

Theorem II. 3 For any one-parameter model, the 
Holevo bound coincides with the SLD CR bound, i.e., 

Cf = -, (30) 

99 

holds for all 9 € 0 where gg is the SLD Fisher informa¬ 
tion at 6. 
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Note that there is no weight matrix since we are dealing 
with a scalar MSE for the one-parameter case. Impor¬ 
tantly, there is no gain from collective POVMs for one- 
parameter models. Existence of a POVM whose MSE is 
equal to this bound is discussed independently by several 
authors [22-24]. 

2. D-invariant model 

Consider an arbitrary /c-parameter model A4 and let 
Lg p (i = 1,2,..., k) be the SLD operators at 8. The lin¬ 
ear span of SLD operators with real coefficients is called 
the SLD tangent space of the model at 8: 

Tq(M) '.= span R {Lg > i, Lg )2 ,..., Lg, k }. (31) 

Any elements of the SLD tangent space, X £ Tg(Xi), 
satisfy tr (pgX) = 0 and it is not difficult to see that the 
space Tg(JVl) is essentially a real vector space with the 
dimension k. Holevo introduced a super-operator V p , 
called a commutation operator, as follows. Given a state 
p on H, let C(T~L) be the set of linear operators on H, 
then Vp is a map from C((H) to itself defined through the 
following equation: 

(Y,V P (X)) P = [Y, X)p, VX,E £ C(H). (32) 

Here, (X,Y) P = Retr(pX*E) is the SLD inner prod¬ 
uct and [A, Y] p = tr {p[X*, E]) /(2i) is a sesqui-linear 
form. (Here, the definition is different from the original 
one by a factor.) When considering a parametric model, 
we denote V pe = Vg for simplicity. We say that a model 
is D-invariant at 8 if the SLD tangent space at 8 is in¬ 
variant under the action of the commutation operator. 
Mathematically, this definition is expressed as 

M is D-invariant at 9 Vg(X) G Tg(M), VX £ Tg(Ai). 

When a model is D-invariant for all 8 £ 0, we say 
the model Xi is globally D-invariant. Lemma B.l in 
Appendix B 2 characterizes equivalent conditions for D- 
invariant models. 

From the definitions of two inner products (9) and the 
commutation operator, the relationship 

(X,Y)+ = (X, {I + W p ){Y))p 

& {X, V p (Y))p = -i«X, E>+ - (A, Y) p ), (33) 

holds for all linear operators A, Y on H. For a given 
model, another important relation 

tr (dipgX) = (A, Lg A )p e = (A, L e ,i)+ , (34) 

holds for VA £ C{H). Combining them gives 

(A ,Lg t i) pe = (A, (/ + iVg)(Lg ^)} pg , and hence we obtain 

L e ,i = {I + iV g )(L e J. (35) 


Two more useful relations are 

(L l g,Vg(L 3 0 )) pe =lmz l e 3 , (36) 

(37) 

which can be checked directly from the definitions. 

It is well known that the Holevo bound gets simplified 
significantly if the model is D-invariant [15]. Importantly, 
D-invariant model enjoys the following proposition, which 
is due to Holevo: 

Proposition II.4 Let Lg (Lg) be the SLD (RLD) dual 
operator and Gg (Gg) be the SLD (RLD) Fisher infor¬ 
mation matrix, respectively. Define a k x k hermite ma¬ 
trix by Zg = [(Dg, Lg)+ ]j,je{i. k}- When the model is 

D-invariant at 9, Zg = G(f l holds at 9 and further the 
Holevo bound is expressed as 

VIE G W Cf [IE] = Cf [IE] = hg\L\W). (38) 

This statement can be proven in several different manners 
[10, 13, 15]. 

In passing, we note that the expression hg[L\W] in 
the above proposition is also expressed as hg[L\W] = 
Tr {WRe Zg} + TrAbs {IEIm Zg} in terms of the matrix 
Zg. When the model is not D-invariant, he[L|IE] does 
not seem to play any important role. This is because the 
quantity 

Cf [IE] := Tr {IERe Zg} + TrAbs {IEIm Zg} , (39) 

is always greater or equal to the Holevo bound, i.e., 
Cf [IE] > C(f [IE] for all weight matrices. Nevertheless, 
as will be shown in this paper, this is an important quan¬ 
tity and we call it as the D-invariant bound in our discus¬ 
sion. We note that this quantity (39) was also named as 
the generalized RLD CR bound by Fujiwara and Nagaoka 
[25] in the following sense. When a model fails to satisfy 
some of the regularity conditions, the RLD operators do 
not exist always. Even in this case, when the model is 
D-invariant, then the above bound (39) is well defined 
and provides the achievable bound for a certain class of 
models, known as the coherent model [25]. 

Another remark regarding this proposition is that the 
converse statement also holds. 

Theorem II.5 For any k-parameter model M. on any 
dimensional Hilbert space under the regularity conditions, 
the following equivalence holds: 

M. is D-invariant at 8. 

VIE G W C^[IE] = Cf [IE], 

VIE G W Cf [IE] = Cf [IE]. (40) 

This equivalence for the D-invariant model might have 
been known for some experts, but it was not stated ex¬ 
plicitly in literature to our knowledge [26]. Sketch of 
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proof is given in Appendix B 3 for the sake of reader’s 
convenience. 

We remark that the Holevo bound for a general model, 
which is not D-invariant, exhibits a gap among Cf[W\, 
and Cg\W}. The following relation holds in gen¬ 
eral: 


c$ [W] > Cf [W] > max{C# [W ], Cf [W] }, (41) 

for all weight matrices W. From this general inequality, it 
is clear that the condition of D-invariance shrinks the gap 
between Cg[W] and Cg-\W\ to zero. The Holevo bound 
then coincides with the RLD and D-invariant bounds. 


III. THE HOLEVO BOUND FOR QUBIT 
MODELS 


In this section we consider a model for quantum two- 
level system, a qubit model. For mixed-state models, 
where parametric states are rank-2 for all 9 € 0, possi¬ 
ble numbers of the parameters are from one to three. As 
stated in Sec. IIC, the Holevo bound for one-parameter 
qubit model is solved and is given by Theorem II.3. When 
the number of parameters is equal to three, on the other 
hand, it is easy to show that the model becomes D- 
invariant as follows. Since three SLD operators are lin¬ 
early independent, they expand the set of all linear oper¬ 
ators X satisfying the condition tr ( pgX ) = 0. In other 
words, the SLD tangent space is same as this space and 
hence the SLD tangent space is always D-invariant. In 
this case, the Holevo bound is given as Eq. (29). There¬ 
fore, the two-parameter case needs to be solved explicitly. 
In the following, we consider two-parameter qubit models 
of mixed states only. Further, the regularity conditions 
mentioned before are assumed throughout our discussion. 

Upon performing this optimization to derive an explicit 
formula for the Holevo bound, it is convenient to utilize 
the Bloch-vector formalism. A similar technique has been 
used by Watanabe, where all operators are expanded in 
terms of a basis of Lie algebras [27]. In the next subsec¬ 
tion, we present necessary machinery and then solve the 
two-parameter case. 


A. Bloch-vector formalism for qubit estimation 
problem 

In this subsection, we present a formalism in which 
SLD operators are represented by three-dimensional real 
vector. This is motivated by the well-known one-to- 
one mapping between a given qubit state and three- 
dimensional real vector. Thus, any qubit model can be 
represented by a family of three-dimensional real vectors 
as 


with B = {x € K 3 | |x| < 1} the interior of the Bloch ball. 
To simplify notations, we define the standard inner prod¬ 
uct and the outer product for three-dimensional complex 
vectors by 

(a\b) = ^2 |a)(b| = 

2 —1,2,3 

respectively, where a denotes the complex conjugation of 
a. The outer product is a 3 x 3 matrix whose action onto 
a vector c £ C 3 is |a)(6|c = (b|c)a. 

We first observe that the one-to-one correspondence 
between the SLD operator and a four-dimensional vector 
when is expanded in terms of the basis {/, <ti, 02 , ( 73 } with 
<jj usual Pauli spin matrices for spin-1/2 particles. Since 
the SLD operators belong to the SLD tangent space, the 
relation tr {pgLgi) = (I,Lg^) pe = 0 holds, i.e., they are 
orthogonal to the identity operator with respect to the 
SLD inner product. This leads to the following con¬ 
straint: 



tr (Lg ti ) = ~{sg\lgg), (43) 


where ig p = l 2 ei , l 3 e i ) T with ^e,i = tr ( a j L o,i) is a 
three-dimensional real vector. Thus, we have a one-to- 
one mapping from the SLD operator Lg p to the three- 
dimensional real vector £g p . The vector £g p shall be re¬ 
ferred to as the SLD Bloch vector in this paper. 

It is straightforward to solve the operator equation 
(10), which defines the SLD operators, and the SLD 
Bloch vector is 


diSg + 


(<9 j; s e |s0) 

1 - s 2 
1 Sg 


s$, 


(44) 


where sg = |sg| denotes the length of the Bloch vector 
and d,j = d/dO 1 is the ith partial derivative. To proceed 
further, we find it convenient to introduce a 3 x 3 matrix: 


:= 1 


|se)(s0| 

1 _ s -2 ' 

± 


(45) 


with 1 the identity matrix acting on the three- 
dimensional vector space C 3 . It follows from the defini¬ 
tion that Qg is a real and positive matrix with eigenvalues 
1 ,1, (1 — s 2 e )~ l and its inverse is 


Q? = l-\a 0 ){a 9 \. (46) 

The SLD Bloch vector is then expressed as 

£g,i = Qgdi.sg (<t=> Q g x igg = diSg). (47) 


The (i. j) component of the SLD Fisher information is 
9 e,ij = = {dis e \Qgdjs e ). (48) 


Let g l g J = ( Gg 1 )y be the (i,j) component of the inverse 
SLD Fisher matrix and we define the SLD Bloch dual 
vector £g by 


ft = 






J 


■Mb — { s e — (sg, Sg, Sg) £ B | 6 £ 0} , (42) 


(49) 



B. Two-parameter qubit model 


then, the following orthogonality condition holds: 

= %> ( 50 ) 

which corresponds to Eq. (13). The inverse of SLD Fisher 
information matrix is also expressed as 

gl j = ( 51 ) 

The same line of arguments holds for RLD operators 

and RLD Fisher information. The only difference is here 
is that the RLD Bloch vector becomes complex in general. 
Define a 3 x 3 complex matrix: 

Qe '■= -——(1 — iFg) , (52) 

1 — s e 

where ( Fg)ij := eikjSg t k with eikj the completely anti¬ 
symmetric tensor. The action of Fg is to give the exterior 
product of two vectors, i.e., Fga = sg x a for a £ C 3 . 
From this definition, Qg is also strictly positive and its 


inverse is given by 

Qg 1 = t-\a e )(a e \+iF e = Qg 1 +iF e . (53) 

The RLD Bloch vector is 

h,i = QediSg (-o- Qg X lg,i = d z sg), (54) 

and the RLD Fisher information matrix is 

9e,ij = (h,i\Qg'&ej) = (' diSe\QgdjSg ). (55) 

Define the RLD Bloch dual vector by 

^e = Y^9e^o,j, (56) 

3 

then we have 

(4l Qg 1 te, j )=S i j . (57) 


Other useful relations are listed below without detail 
calculations. First, there is a one-to-one correspondence 
between SLD and RLD Bloch vector. This is given by 

ie,i = (1 + i Fg)£g,i (■» ie ti = (Qg 1 - i F e ) tg ti ). (58) 

Second, the vector £g,i~£g,i is orthogonal to the Bloch 
vector sg. Defining 

le,i ■= (se\t-e,i), le,i ■= (sg\£g,i), (59) 

this is expressed as 

le,i = 18,i- (60) 

Third, the SLD Fisher information and real part of the 
RLD Fisher information are related by 

2 98,ij ~ Re gg tij = i8,ne,j■ (61) 

1 — s e 

In other words, the matrix (1 — s^) _1 Gg — Re Gg is rank 
one. 


In this subsection, we consider an arbitrary two- 
parameter qubit model, that is the parameter to be es¬ 
timated is 9 = (0i,0 2 ) 6 0. In order to derive an ex¬ 
plicit expression for the Holevo bound, we first rewrite 
the Holevo bound (24) in terms of the Bloch vectors. 

A linear operator which satisfies tr ( pgX *) = 0 can be 
expressed as 

X i = -(sg\x i )I + X i -(T. (62) 

Let Tgh = {x £ R 3 | (x\d-iSg) = 0} be the orthogonal 
space to the ith derivative of the Bloch vector and an 
element of the set appearing in the definition (20) takes 
the form of A = (A 1 , A 2 ) with 

A* = —(sg\x l )I T x l • cr , 

x l £ Tg-j for j ^ i and (x l \diSg) = 1. (63) 

Thus, the set of operators X can be mapped to the set 
of vectors x = ( x 1 ,x 2 ) T £ R 6 . Using this form of Bloch 
vector representation, the (i,j) component of the Zg[x\ 
matrix and the Holevo function read 

Zg[x\ = (x^Qg'-x 3 ), (64) 

2 

hg[x\W} = ^ | w i j(x l, \Qg 1 x-’) + Vdet W\(x l \Fgx-'}\ , 

i,j— i 

for a given 2x2 weight matrix W = [iUy]i,je{i, 2 }- 

We note that the Holevo function (64) is a quadratic 
function of the six-dimensional vector x. The minimiza¬ 
tion of this function under the constraints (63) can be 
handled with a standard procedure. The only point 
needs to be taken is that the function is not differen¬ 
tiable for all points. Since the number of free variables 
for the optimization is 6 — 4(the number of independent 
constraints)= 2, we take the following substitution: 

** = 4+(65) 

where ij- = d±sg x <9 2 sg is a vector orthogonal to both 
disg and 9 2 se and £ = (f 1 ,f 2 ) T £ R 2 is a free vari¬ 
able without any constraint. With this substitution, the 
Holevo function is significantly simplified as follows. 

Lemma III.l For a two-parameter qubit model, the 
Holevo bound takes the following minimization form with¬ 
out any constraint: 

C?[W] = min hg[£\W], ( 66 ) 

{SR 2 

where the function /ie[£|W] is defined by 

hg [Cl W] = Tr {WGg 1 } + (^IQ^Hljuf) 

+ 2v / detTU (iWFgl 2 ) + (1 - s 2 e )mo\ . (67) 
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In this expression, we introduce the standard inner prod¬ 
uct for two-dimensional real vector space by (a\b) = 
aih + a 2 b 2 and 70 = (jg, 1 , Tg, 2 ) T is given by Eq. (59). 
The derivation for this lemma is given in Appendix A 1. 


C. Main result 

In the following, we carry out the above optimization 
to derive the main result of this paper, Theorem 1.1. We 
first list several definitions and lemmas. 

Definition For a given two-parameter qubit mixed-state 
model, the SLD CR, RLD CR, and D-invariant bounds 
are defined by 

C$[W]=Tv{WG^}, (68) 

C$[W\ = TV {iTReG'-'} + TrAbs {VlmG ^ 1 j , 

C$\W] = Tr {WRe Zg} + TrAbs {WTm Zg} , 

respectively, where Gg is the SLD Fisher information 
matrix, Gg is the RLD Fisher information matrix, and 

Zg ■= [ z e]i, je{ 1 , 2 } with z e : = tr ( Pe L g L g ) as before. 

Lemma III.2 For any two-parameter qubit mixed-state 
model, the following relations hold: 

1. (£g\ = (1 - s|))det Gg = (1 - si) 2 det Gg. 

2. 2Vdet W | (£l\Fg£g) \ = TrAbs jlLTmG^ 1 } = 

TrAbs {WIm Zg}. 

3. mw~%) = de \ W \ Ge (C$[W] - C?[W}). 

1 Sg 


(0,0,0), the model is always D-invariant at this point. 
This is because the condition 7 g ti = ( sg\£g ti ) = 0 is met 
at sg = (0,0,0). Last, a globally D-invariant model is 
possible if and only if the state is generated by some uni¬ 
tary transformation. This is because the condition III.3-6 
in Lemma III.3 is equivalent to preservation of the length 
of the Bloch vector. 

Finally, we need the following lemma for the optimiza¬ 
tion: 


Lemma III.4 For a given positive 2x2 matrix A, a real 
vector b £ R 2 , and a real number c, the minimum of the 
function 


m = (t\M )+2 m)+ c 


is given by 


min /(I) 

££R 2 


2|c| — (b\A 1 b) if |c| > (b\A 1 b) 


(b\A~ l b) 


if |c| < {b\A- lz b) 


where the minimum is attained by 

f — sign(c)A _1 6 if |c| > ( 6 |A _ 1 6 ) 

| —-— —A^b if |c| < ( 6 |A _ 1 6 ) 

{ (b\A~'b) 

where sign(c) is the sign of c. 

Proofs for the above three lemmas are given in Appendix 

A2-A4. 


1. Proof for Theorem 1.1 


Lemma III.3 For any two-parameter qubit model A4, 
the following conditions are equivalent. 

1. JA is D-invariant at 9. 


We now prove Theorem 1.1. From the expression of 
the Holevo function (67), we can apply Lemma III.4 by 
identifing 


2. ReG 0 1 = G g 1 at 9. 

3- 70,1 = 70,2 = 0 at 9. 

Furthermore, we have the following equivalent character¬ 
ization for global D-invariance. 

f. M. is globally D-invariant. 


A = (£g~\Qg 1 £g~)W, 

(69) 

b= (1 - Sg)YdetWjg, 

(70) 

c = YdetW(£ 1 e \Fg£l). 

(71) 


We need to evaluate (b\A 1 b) and |c| and they are calcu¬ 
lated as follows. 


5. ReG^ 1 = G g 1 for all 9 £ 0. 

6 . |s@| is independent of 9. 

Three remarks regarding the above lemmas are listed: 
First, imaginary parts of the inverse of the RLD Fisher 
information matrix and the Zg matrix are always iden¬ 
tical for two-parameter qubit mixed-state models, i.e., 
ImG ^" 1 = Im Zg, see proof in the appendix. Second, if a 
model expressed as the Bloch vector contains the origin 


(b\A~ l b) = (l-s1YdetWm((K\Qgtg)W)-^g) 


(1 — Sg) 2 det W 


(■lg\Qg X £g) 

= C Z g\W) - C?[W] 


(le\W 1 Te) 


where Lemma III.2-1 and III.2-3 are used to get the last 
line. Lemma III.2-2 immediately gives 

2|c| = TrAbs {WlmZgj = C$[W] - Cjj[W\. (72) 
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Therefore, we obtain if 

|c| > (b\A-'b) & C?[W] > Uc$[W] + Cf\W]) 


which is continuous and the first derivative is also con¬ 
tinuous for all x £ R. Then, the Holevo bound is written 
in a unified manner: 


is satisfied, the Holevo bound is C H \W] 

Cg[W] = C$[W} + {Cg [IT] - Ci[W\) - {Cg [IT] - C?[W\) 

= Cf[W]. (73) 


C$\W\ + {Cf[W\ - C?[W]) 


x H 


(C$[W) - C$[W})/2 \ 

c§\w\-c*[W] ) 


(79) 


If |c| < ( 6 |A _ 1 6 ) 47 Cg[W\ < (Cg[W] + Cg\W])/2 is 
satisfied, on the other hand, the Holevo bound takes the 
following form: 


C?[W] 


Cg[W] + 
C?[W] 


[(C?[W]-C§[W])/2] 2 
C z e [W} - C*[W\ 


i i ( c$m + cjm 
C?[W}-C*m V 2 
= Cf[W] + Sg[W], 



where the function 5e[W] is defined in Eq. ( 6 ). This 
proves the theorem. □ 

We remark that from Lemma III.2 and the positivity 
of IT we always have the relation 

C?[W]>C?[W], (75) 

and the equality if and only if {'fe\W~ 1 ^g) = 0 47 M 
is D-invariant at 9 by Lemma III.3-3. Note if M. is D- 
invariant at 9 {"ie,i = 7 g i2 = 0), the condition [IT] < 
(Cf[W\ + C$[W\)/2 « 0 < TrAbs{ITIm Gg 1 }/2 < 
Cf[W] — (^[ir] cannot be satisfied. Thus, the obtained 
Holevo is well defined for all 9 and for arbitrary weight 
matrix IT. 

The optimal set of hermite operators attaining the 
Holevo bound can be given by Lemma III.4 as follows. 
Define an hermite matrix by 

L e '■= -( s e\te) I + te ' ( 76 ) 


and the function £(IT) by 


(e{W) := AetW-'^Gg 1 
I sign(Im Zg 2 ) if C?[W\ > \{C Z [W] + C$[W}) 

X \ W^Imzh 2 

[cfWTcfW] otherwise 
Then, we have X * = {Xf,X 2 ) = axgmm.hg[X\W\ as 
Xl = Lg + <^g{w22le,i — 1*7270,2 )L'o > 

X 2 = Lg + Ce(~W 2 i"fe,i + 1 * 7170 , 2 )Lg. (77) 

Before moving to discussion and consequence of the 
main result, we present the following two-alternative ex¬ 
pressions. Define a function, 


H{ x) 


2 \x\ — 1 if jar| > 1 
x 2 if \x\ < 1 


(78) 


This expression needs a special care when C Z [W] — 
Cg- [IT] —> 0. This case should be understood as the limit 
lim^o aH(b/a) = 2|6|. 

The other expression shown in Eq. (7) follows 
from the first line of Eq. (74) by noting C Z [W] — 
Cf [IT] = Tv{W{G^ -ReGg L )} and C$ [IT] - C§ [IT] = 
Tr Abs{ ITIm Gg 1 }. 


IV. DISCUSSION ON THEOREM 1.1 

In this section, we shall discuss the consequences of 
Theorem 1.1. This brings several important findings of 
our paper. First is two conditions that characterize spe¬ 
cial classes of qubit models. Second is a transition in the 
structure of the Holevo bound depending on the choice 
of the weight matrix. 


A. Necessary and sufficient conditions for special 
cases 

The general formula for the Holevo bound for any two- 
parameter model is rather unexpected in the following 
sense. First of all, it is expressed solely in terms of the 
three known bounds and a given weight matrix. Sec¬ 
ond, a straightforward optimization for a nontrivial func¬ 
tion reads to the exactly same expression as the RLD CR 
bound when the condition Cg' [IT] > (Cg\W] + Cg\W])/2 
is satisfied. As noted before, this condition explicitly de¬ 
pends on the choice of the weight matrix IT. At first 
sight, this seems to be in contradiction with the general 
theorem II.5 stating that the RLD CR bound can be at¬ 
tained if and only if the model is D-invariant. Therefore, 
we must examine that the Holevo bound is identical to the 
RLD CR bound if and only if the model is D-invariant 
based on Theorem 1.1. The following proposition con¬ 
firms that this is indeed so. We note that this statement 
is a special case of Theorem II.5. Here, its proof becomes 
extremely simple with the obtained formula. 

Proposition IV. 1 For any two-parameter qubit model, 
the Holevo bound at 9 becomes same as the RLD CR 
bound for all positive weight matrices if and only if the 
model is D-invariant at 9. That is, 

VIT £ W Cg?[W\ = Cg' [IT] 47 M is D-invariant at 9. 

(80) 
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Proof: When the model is D-invariant, the condition 
Cq-\W\ < (Cf[W] + Cg[W])/2 cannot be satisfied. (See 
the remark after Eq. (75).) Therefore, the Holevo bound 
is always identical to the RLD CR bound in this case. 

Next we show the left condition implies the right in 
Eq. (80). If Cg 1 [IE] = C/* [IE] holds for all IE, the relation 
VdetW | (£l\Fg(.g)\ > Tr{IE(G^' 1 — ReG/ -1 )} must be 
true for all IE. But, if Gg 1 — Re (Gg 1 ) ^ 0, this is impos¬ 
sible. To show it, let us suppose B := Gg 1 — ReG/" 1 > 0, 
then we can change IE — > B~ 1 ^ 2 W'B~ 1 ^ to get 


Vdet W , i „ „2\ 


Tr (IE'} 


{tg\Fgtg)\ det B- 1 > 1, VIE' G W. (81) 


It is easy to show that /(IE') = \J det lE'/Tr {IE'} as a 
function positive matrix IE' satisfies 1/2 > /(IE') > 0 as 
changing positive matrix IE'. That is /(IE') can be made 
arbitrary small by choosing IE' and hence the condition 
(81) cannot hold. This gives B = 0 <G> Gg 1 = ReG^ 1 
and hence this is equivalent to the D-invariant condition 
from Lemma III.3. □ 


Using Theorem 1.1, we now state one more impor¬ 
tant condition characterizing the model where the Holevo 
bound coincides with the SLD CR bound. 

Proposition IV.2 For any two-parameter qubit model, 
the Holevo bound coincides with the SLD CR bound for 
all positive weight matrices if and only if the imaginary 
part of Zg vanishes. Further, this condition is equiva¬ 
lent to the existence of a weight matrix Wq such that the 
Holevo bound and the SLD CR bound are same with this 
particular choice IEq [30]. Mathematically, we have 


VIE G W Cf [IE] = Cg\W\ (82) 

<G> Im Zg = 0 (83) 

3 Wo G W Cg [Wo] = Cg [Wo] . (84) 

Proof: (83) => (82): The condition Im Zg = 0 <G> 


Im (Lg^i, Lg^ 2 )+ g = 0 implies TrAbs{WImG e 1 } = 0 for 
all W G W. In this case, the condition G<f[IE ] < 
{Cg\W\+ Cg\W])/2 is always satisfied, and hence the 
Holevo bound coincides with the SLD CR bound for all 
choices of the weight matrix. This follows from the sec¬ 
ond line of the expression (74). 

(82) => (84): The condition (82) implies an existence 
of a weight matrix satisfying G^ [IE] = Cg [IE ]. 

(84) => (83): Next let us assume C^flEo] = 
C$\W 0 ] for some weight matrix Wo G W. The ex¬ 
pression (79) immediately implies TrAbs {Wolrn Zg) = 
2-^/det Wo|Im (Lg t i, Lg )2 )+ | = 0 which gives Im Zg = 0. 
Therefore, three conditions are all equivalent. □ 

We have three remarks regarding this proposition. 
First, in terms of the SLD Bloch vectors, the necessary 


and sufficient condition (83) is also written as 
Im Zg = 0 <G> tr {pg [Lgi, Lg , 2 ]) = 0 
<=> (sg\Lg t i X £g^) = 0 
<G> {sg\disg X d 2 sg) = 0, 

which is easy to check by calculating the Bloch vector of 
a given model. 

Second, we note that given a symmetric matrix A, 
Tr{WA} > 0 for all positive weight matrix W implies 
A > 0 as a matrix inequality. When the Holevo bound 
is same as the SLD CR bound, we see that the MSE 
matrix for any asymptotically unbiased esti¬ 

mators satisfy the SLD CR inequality: 

lim nVg {n) [fl^\ > Gg 1 , (85) 

n—>• oo 

as a matrix inequality. Moreover, there exists a sequence 
of estimators that attains this matrix equality. This is 
rather counter intuitive since two SLD operators Lg i and 
Lg 2 do not commute in general. Therefore, the condition 
Im (Lg t i, Lg t2 )+ g = 0 seems to grasp asymptotic commu¬ 
tativity of two SLD operators in some sense. Indeed, this 
condition can be written as tr (pg [Le.i, Lg ^]) = 0, i.e., 
commutativity of Lg i and Lg 2 on the trace of the state 
pg. When this holds, the quantum parameter estimation 
problem becomes similar to the classical case asymptoti¬ 
cally. In the rest of the paper, we call a model asymptot¬ 
ically classical when this condition is satisfied. A similar 
terminology, “quasi classical model,” was used by Mat- 
sumoto in the discussion of parameter estimation of pure 
states [28]. Here, we emphasize that classicality arises 
only in the asymptotic limit and hence, this terminology 
is more appropriate. We also note that the equivalence 
between (83) => (84) was stated in the footnote of the 
paper [9] based on the unpublished work of Hayashi and 
Matsumoto. Here our proof is shown to be simple owing 
to the general formula obtained in this paper. 

Last, a great reduction occurs in the structure of the 
fundamental precision bound for this class of models. We 
note that achievability of the SLD CR bound for specific 
models have been reported in literature in discussion on 
quantum metrology [31, 32]. Here, we provide a simple 
characterization, the necessary and sufficient condition, 
of such special models in the unified manner. 

Having established the above two propositions, we can 
conclude that a generic two-parameter qubit model other 
than D-invariant or asymptotically classical ones exhibits 
the nontrivial structure for the Holevo bound in the fol¬ 
lowing sense: The structure changes smoothly as the 
weight matrix IE varies. For a certain choice of IE, it co¬ 
incides with the RLD CR bound and it becomes different 
expression for other choices. Put it differently, consider 
any model that is not asymptotically classical, then we 
can always find a certain region of the weight matrix set 
W such that the Holevo bound is same as the RLD CR 
bound. This point is examined in detail in the next sub¬ 
section and examples are provided in the next section for 
illustration. 
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B. Smooth transition in the Holevo bound 

Let us consider a two-parameter qubit model that 
is neither D-invariant nor asymptotically classical. In 
this case, the set of all possible weight matrices 
W = {W\W is 2 x 2 real positive definite} is divided 
into three subsets. The first two sets are W + and W_ 
in which C<f[W] — ( Cg\W} + Cg\W])/2 is positive and 
negative, respectively. The last is the boundary Wg that 
consists of a family of weight matrices satisfying the equa¬ 
tion: 

B e [W] := C*[W] - \(C z e \W] + C§[W]) = 0. ( 86 ) 

According to this division, the Holevo bound takes the 
form C^[W] = [IF] for W G W+ and it is expressed 
as C^[W] = C$[W\ + S e [W] for W G VV_. Whereas, 
C(f\W] = Cg^IF] = Se[W] holds for the boundary W G 
Wg. 

In the following we characterize these sets explicitly. To 
do this, we first note that the degree of freedom for the 
weight matrix in our problem is three due to the condition 
of IF being real symmetric. Second, we can show that a 
scalar multiplication of the weight matrix does not change 
anything except for the multiplication of over all expres¬ 
sion of the Holevo bound. Thus, we can parametrize the 
2x2 weight matrix by two real parameters. For our pur¬ 
pose, we employ the following representation up to an 
arbitrary multiplication factor: 


where the common conditions w 2 > 0 , |io| < 1 , and c > 0 
also need to be imposed to satisfy the positivity of the 
weight matrix. 

C. Limit for pure-state models 

So far we only consider models which consist of mixed 
states. It is known that collective measurements do not 
improve the MSE for pure-state models [4, 13, 28]. In 
other words, the Holevo bound is same as the bound 
achieved by separable measurements as far as pure-state 
models are concerned. 

In this subsection, we examine the pure-state limit for 
our general result. When a mixed-state model is asymp¬ 
totically classical, the Holevo bound is identical to the 
SLD CR bound. This should be true in the pure-state 
limit, and this agrees with the result of Matsunroto [28]. 
When a model is D-invariant, on the other hand, it is 
shown that the RLD CR bound can be achieved. This 
also holds in the pure-state limit and we examine the 
pure-state limit for a generic mixed-state model below. 

We first note that we cannot take the limit for SLD and 
RLD operators directly [25, 29]. This is because there are 
the terms 1 — Sg appearing in the denominators. However, 
the SLD and RLD dual operators are well defined even in 
the pure-state limit. By direct calculation, we can show 
that the SLD and RLD dual vectors (49, 56) are written 
as 


IF = Ug 


1 WW 2 
WW 2 w\ 


US 


e i 


(87) 


where w 2 > 0 and detIF = w 2 (l — w 2 ) > 0 => \w\ < 1 
are imposed from the positivity condition and the real 
orthogonal matrix Ug is defined in terms of Eq. (59) by 


Ug 



76.1 — 76,2 

76.2 76,1 


( 88 ) 


The assumption of the model under consideration then 
yields 


Bg[W] = \{l l g \Fgi 2 g )\VAetW — Tr |lT(G ^ 1 — ReG^ 1 ) j 


\{l\\F e tl)Wl - 1 v 2 ~ Tr {G ^ 1 - ReGg 1 j w 2 


W 2 - 


Therefore, by defining a constant solely calculated from 
the given model ag := K^l-Fe^l/Tr jG )]" 1 — ReG^ 1 !, 
we obtain the sets W±, W 3 as follows. 



1 agww 2 

agww 2 o? s w\ 

1 agww 2 

agww 2 OgW 2 





1 agww 2 
agww 2 agW 2 


U*g 

U* e 

Ug 


W 2 + W 2 < 1 } , 


w\ + W 2 > 1 } , 


wl + w 2 = 1 } , 


(89) 


Qg X lg X d 2 Sg 

<#l QeX) ’ 

Q~g X tg X dl Sg 

QoX) ’ 

(1 - i Fg)(£j- x d 2 sg - i(se|lg- )d 2 sg) 
(l^Qg^g) 

(1 - \Fg)(ljf X di Sg - \{sg\l^)d!Sg) 

{£g\Qg lf -g) 

for any two-parameter qubit model. Thus, as long as 
(tj-\Qg £g) 7 ^ 0 , they converge in the pure-state limit. 
This condition is also expressed as \ljf x s @| 2 — (1 — 
s ^) 2 |£^| 2 ^ 0 , and hence is equivalent to if x sg ^ 0 in 
the pure-state limit., i.e., M is not asymptotically clas¬ 
sical. Second, the same warning applies to the SLD and 
RLD Fisher information [25, 29]. However, the inverse of 
the SLD Fisher information matrix is well defined even 
for the pure-state limit. This is because the (i,j) compo¬ 
nent of the inverse of SLD Fisher information matrix is 
{f-glQ^ 1 Ug) and this has the well-defined limit. The same 
reasoning can be applied to the inverse of the RLD Fisher 
information matrix in the pure-state limit. 

Last, let us examine the general formula. It is straight¬ 
forward to show that the function ( 6 ) vanishes in the 
pure-state limit. In other words, the general formula 
given in Theorem 1.1 becomes C^[W] = Cg^W] for all 
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weight matrices. Therefore, the Holevo bound becomes 
the RLD CR bound in the pure-state limit. This is, of 
course, expected because any pure-state model can be 
expressed as a unitary model, which is D-invariant from 
Lemma III.3. 


V. EXAMPLES 


In this section, we consider several examples for two- 
parameter qubit models to illustrate our result. The first 
one is a D-invariant model whose Holevo bound is iden¬ 
tical to the RLD CR bound. The second one is a asymp¬ 
totically classical model which gives the SLD CR bound. 
As the last example, we analyze a generic model in par¬ 
ticular the behavior of the Holevo bound depending on 
the weight matrix. 

Within the setting of pointwise estimation, we note 
that it is sufficient to specify the following three vectors 
to define a model locally; 

{so, disg, d 2 sg} at0. (90) 

Equivalently, the set {sg,£gp,£g 2 } or {.sg, £e,i, £ 9 , 2 } can 
be used to define the model at 6 uniquely. One-to-one 
correspondences among these three specifications of the 
model are easily established by the Bloch vector repre¬ 
sentation discussed in Sec. Ill A. 

The D-invariant condition is now expressed in terms of 
the vectors of the set (90) as 

D-invariant <t=> (s^c^s#) = (sq^sq) = 0. (91) 


B. Asymptotically classical model 


Consider the following model: 

Mb = {a g = f 1 (0)t*i + f\e)u 2 | 6 e 0} , (94) 

where it; are unit vectors, which are not necessarily or¬ 
thogonal to each other, and / J '(0) are scalar (differen- 
tialble) functions of 9 = (0 3 ,0 2 ). The parameter region 
0 is specified by an arbitrary open subset of the set; 
{(0i, 62 ) | sg < 1}. We can show that this model is asymp¬ 
totically classical because of x d 2 sg oc u\ x u 2 , 

and the Holevo bound is C^\W] = Tr {WGg 1 } MW > 0. 

For the asymptotically classical case, we can easily 
compare this result with the bound achieved by an op¬ 
timal estimator comprised of separable POVMs. In this 
case, the Nagaoka bound is known to be achievable that 
is calculated as [21] 


Cg [W] = Tr {WGg 1 } + 2 sJ&etWGg 1 . (95) 


Therefore, the gain due to use of collective POVMs is 
2y / det WGg 1 . 

As a special case, let us set u 3 to be orthogonal normal 
vectors and P( 6 ) = 0A Then, the inverse of the SLD 
Fisher information matrix reads 


1 _(l- ( 0 1 ) 2 - 9 3 0 2 \ 

W ^ -0 l 9 2 1 - (0 2 ) 2 ) ’ 


(96) 


and the gain mentioned above is 2y/det W\J 1 — Sg. Thus, 
the role of collective POVMS becomes important as the 
state becomes more mixed. 


The asymptotically classical condition is equivalent to 

Asymptotically classical <t=> (s^Sis# x d 2 sg) = 0. (92) 

All models other than satisfying the above two conditions 
are generic ones. 


A. D-invariant model 

Consider a simple unitary model where the parametric 
state is generated by a two-parameter unitary U{9) (9 = 
( 01 , 02 )): 

pg = U(9)p 0 U(9)*. (93) 

It is easy to show that the norm of the Bloch vector is 
independent of 9 , and hence, this model is globally D- 
invariant. The Holevo bound is Cg [IT] = C^\W] for all 
weight matrices W. 

We note that this kind of D-invariant models possesses 
symmetry and has been studied by many authors, see for 
example Ch. 4 of Holevo’s book [15] and Hayashi [7] and 
references therein. 


C. Generic model 

In this subsection, we analyze a generic model other 
than the previous two examples. In this case the struc¬ 
ture of the Holevo bound changes when the weight matrix 
W varies, that is it takes the same form as the RLD CR 
bound for certain choices of W whereas it becomes differ¬ 
ent form for other cases. This is explicitly demonstrated 
in the section IV B and here we consider a simple yet 
nontrivial example to gain a deeper insight into this a 
phenomenon of this transition. 

Consider a two-parameter qubit model: 

M B = {sg = (9\ 9 2 ,9 0 ) T \9ee}, (97) 

where 0o is a fixed parameter which satisfies 0 < 

10o I < 1- The parameter region is specified by a sub¬ 
set of {(0 1 ,0 2 )|(0 1 ) 2 + (0 2 ) 2 < 1 — 9q}. This model 
is neither D-invariant nor asymptotically classical when 
0o 7^ 0 and 9 1 9 2 ^ 0 are satisfied. This model be¬ 
comes asymptotically classical for 0 O = 0 as discussed 
before. We can also regard this model as a sub¬ 
model embedded in the three-parameter qubit model: 
{s e = (0 1 ,0 2 ,0 3 ) T |0G0CK 3 }. 
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The inverse of SLD and RLD Fisher information ma¬ 
trices at 9 are 


G,- 1 

G,- 1 


1 -ol 

1 - s e 

1 -ol 


i — — ( 0 1 ) 2 -e l e 2 
-e 1 e 2 i-el- (6» 2 ) 2 

1 —i9q \ 

i^o i ) 


(98) 


with se = [(0 1 ) 2 + (0 2 ) 2 + 91] 1 / 2 , respectively, and the 
three bounds appearing in Theorem 1.1 read 

C$[W] = Tr {W} - WO), 

c * [w] = t^ 0 Ti ' {w}+ 2 ^r^ 0 ldolV ^ W ’ 

C$[W] = Tr {W} - T ± ¥ J0\W0) + 2^1—| %\VtetW, 


where 1 9) = {0 1 ,0 2 ) T is introduced for convenience. 

In the following let us analyze the structure of the 
Holevo bound using the following parametrization of a 
weight matrix: 


with R(u) = ( C ° SW _SmW ) , (99) 
v ' y smw cos a; J v ' 

where we normalize the trace of W to be one and w £ 
(—1,1), w £ [0, 27t). This parametrization is different 
from the one analyzed in Sec. IV B, yet is convenient 
for the purpose of visualization. It is easy to see that 
the effect of the matrix R(ui) is to mix two parameters 
9 1 and 9 2 by rotating them about an angle w. Since 
det W = (1 — u> 2 )/4 for this particular parametrization, 
we see that the RLD CR bound is independent of the 
weight parameter w. The other bounds depend on two 
parameters ( w,ui ). 

We are interested in how the Holevo bound C^[w,oj] 
changes when we vary the weight parameters w, w. Upon 
plotting, we fix the direction of the estimation parameter 
9 as 9 — \0\(l,l)/y/2 with |0| = (0I#) 1 / 2 . We plot the 
Holevo bound for two sets of the state parameters; (a) 
0 O = 0.2 ,9 = 0.346(1,1)^/2 and (b) 9 0 = 0.275 ,9 = 
0.476(1, l)/\/2 for illustration. Figure 1(a) shows the 
Holevo bound for the state parameter (a) as a function of 
the weight-matrix parameter w,u>. In this plot, the gray 
areas indicate the case for Cjf\W] = C^\W] + Sg\W], 
whereas the white-meshed region indicates the case for 
G^[1F] = Cg-\W]. We also show the other state param¬ 
eter setting (b) in Fig. 1(b) with the same convention. 
In both figures we observe smooth transitions between 
two-different expressions discussed in Sec. IV B. 

From these figures, we see that the Holevo bound coin¬ 
cides with the RLD CR bound for relatively large choice 



FIG. 1. The Holevo bound for the state parameter (a) 6q = 
0.2, 0 = 0.346(1, l)/\/2 and (b) Bo = 0.275, 6 = 0.476(1, l)/\/2 
as a function of the wight-matrix parameter ( w,u ) given in 
Eq. (99). The gray areas indicate the case for C^\W\ = 
C^[1F] + S'g[lF], whereas the white-meshed region indicate 
the case for Cg\W] = Cf[W]. 


of the weight matrix in Fig. 1(a), where as the opposite 
observation holds for Fig. 1(b). To gain a deeper insight, 
let us calculate the value of quantities yey := (se\t-e,i) 
(■i = 1,2). Then, we get 70,1 = 70,2 = 0.292 for 
the case (a) and 70,1 = 70,2 = 0.483 for the case (b). 
Since vanishing of these quantities is equivalent to the 
D-invariance of a model (Lemma III.3), we naively ex¬ 
pect that the smaller values of them imply that a model 
behaves more D-invariant-like. Indeed, the examples pre¬ 
sented here agree with our intuition, yet more detail anal¬ 
ysis is needed to make any conclusion. 


D. Transition in the parameter space 


We briefly discuss another important observation of 
this paper. A generic model, other than special cases 
discussed before, exhibits a transition in the structure 
of the Holevo bound for a fixed weight matrix when we 
change the estimation parameter 9. A rough sketch of 
this argument is that a change in the weight matrix is 
amount to that in the parameter and vice versa. This 
is a well-known fact in the pointwise estimation setting 
[25, 33]. Below we briefly show such an example. The 
model is same as the generic model (97) analyzed the 
previous subsection. 

As noted before, a change in the weight-matrix pa¬ 
rameter wo —> wo + w is equivalent to rotate the pa- 
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rameter (9 1 ,9 2 ) by the angle ui. Depending upon the 
choice of the weight matrix, we can also observe a sim¬ 
ilar transition when we change the parameter 9 in the 
set 0. We set the weight matrix to be a diagonal ma¬ 
trix Wo = diag(0.55,0.45) and 9q = 0.35. Figure 2 plots 
the Holevo bound as a function of the state parameter 
9 = ( 9 1 1 6 2 ). The gray area indicates the region where 
Cff\Wo\ = Cg-\W o] + <Se[Wo] holds, whereas the white- 
meshed region indicates the case for C^\Wo\ = C<f[hFo]- 
This shows that the Holevo bound coincides with the 
RLD CR bound for a certain subset of the parameter 
space. 



FIG. 2. The Holevo bound for the weight matrix Wo = 
diag(0.55, 0.45) and 9o = 0.35 as a function of the state pa¬ 
rameter 9 = ( 9 1 ,9 2 ). The gray area indicates the case for 
Cf\W] = C<f [IF] + Se\W], whereas the white-meshed region 
indicates the case for C$\W] = C^[W]. 


VI. CONCLUDING REMARKS 

In this paper, we have derived a closed formula for the 
fundamental precision bound, the Holevo bound, for any 
two-parameter qubit estimation problem in the pointwise 
estimation setting. This bound is known to be asymptot¬ 
ically achievable by the optimal sequence of estimators 
consisting of jointly performed measurements under some 
regularity conditions. Since there exist explicit formu¬ 
las for the Holevo bound for the pure-state qubit model, 
qubit mixed-state models with one and three parameters, 
our result completes a list of the fundamental precision 
bounds in terms of quantum versions of Fisher informa¬ 
tion, which is calculated from a given quantum paramet¬ 
ric model, as far as qubit models are concerned. The ob¬ 
tained formula shows several new insights into the prop¬ 
erty of the Holevo bound for quantum parameter estima¬ 
tion problems. In the following we shall list concluding 
remarks together with outlook for future works. 

First, the necessary and sufficient conditions for the 
asymptotic achievability of the SLD and RLD CR bounds 
are derived when estimating any two-paranreter family of 
qubit states. In particular, the notation of asymptotically 
classical models is proposed, in which all SLD operators 


commute with each other on the trace of a given para¬ 
metric state. In this case, the weight matrix can be elim¬ 
inated and the problem becomes similar to the classical 
statistics in the asymptotic limit, where the SLD Fisher 
information plays the same role as the Fisher information. 
We note that the notion of asymptotic classicality can be 
extended to any models on any finite-dimensional system 
and the same statement obtained for the qubit case holds: 
The Holevo bound coincides with the SLD CR bound if 
and only if the model is asymptotically classical. The de¬ 
tail of this result shall be given in the subsequent paper. 

Second, the RLD CR bound is shown to be achieved 
for a certain choice of the weight matrix even though 
the model is not D-invariant. This result emphasizes the 
importance of the RLD Fisher information for general 
qubit parameter estimation problems. Since the imagi¬ 
nary parts of the inverse of the RLD Fisher information 
matrix and the Zg matrix are always identical for two- 
parameter qubit mixed-state models, we cannot immedi¬ 
ately conclude that this is the general statement or not. 
It might well be that the only real part of the inverse 
of RLD Fisher information matrix plays an important 
role in higher dimensional systems. This deserves fur¬ 
ther studies and shall be analyzed as an extension of the 
present work. 

Third, our result also provides the (unique) minimizer 
to the optimization problem appearing in the definition 
of the Holevo bound. This set of observables, which are 
locally unbiased in the sense of Eq. (16), can be used to 
construct an optimal sequence of POVMs that attains 
the Holevo bound asymptotically. This line of approach 
has been proposed by Hayashi, who reported a theorem 
without proof to realize the construction of POVMs on 
the tensor product states [6]. His approach is different 
from other approaches given in Refs. [10, 11, 13] and it 
may need more refined arguments to complete his theo¬ 
rem. 

Last, smooth transitions in the structure of the Holevo 
bound is shown to occur in general. Since this happens in 
the simplest quantum system, we expect similar phenom¬ 
ena occur in higher dimensional systems as well. How¬ 
ever, we do not know whether or not the number of differ¬ 
ent forms is always two as demonstrated here. The tech¬ 
niques used in this paper can be applied to two-parameter 
estimation problems in higher dimensional systems, and 
we shall make progress in due course. 
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Appendix A: Proofs for lemmas 


1. Leema III.l 

Proof follows from a straightforward calculation. We 
substitute x l = £ g +^£j- into the second line of Eqs. (64). 
The first term reads 

<*‘l Qe 1 ^) = (e e \Qi 1 4)+??(*9\ Qe^e) 

= 9? + (Al) 

where the relation {£ l g \ Q g £ g ) = 0 for i = 1,2 and 
Eq. (51) are used. Note (a\Fga) = 0 for all a £ R 2 , 
then the second term is calculated as 

2 

E \{x i \Fox*)\=2\{x 1 \F e x*)\ (A2) 

*)j= 1 

= Wlm 2 e) + e(4\Fe4) + 

= ml\F e £l) + ^(1 - sl) le ,i + e 2 (l - si) 70, 2 |, 

where {ljf\Fg£ 2 g ) = (1 - and {£^\Fg£\) = -(1 - 

Sg)^ie ,2 are used. Combining the above expressions, we 
get the expression of this lemma. □ 


2. Lemma III.2 

1 . <#|= (1 - s 2 ) det Gg = ( 1 - s 2 ) 2 det G e : 

For convenience, let us define 

n e := disg x <9 2 S0, (A3) 

mg : = £g, 1 X £ 0 , 2 , 

which are related by ng = (1 — s g )Qgmg. The standard 
vector analysis shows 

( m g \mg ) = (£( 9 , 11 ^ 0 ,i)(^0,2^0,2} — ((^ 9 , 11 ^ 0 ,2)) 2 , 

Sg X m e = 70,2^0,1 — 70,1^0,2- 

From the expression ( 48 ), the components of the 
SLD Fisher information matrix is expressed as gggj = 
(£e,i\£g,j) — 70,?70j- Then, the determinant of the SLD 
Fisher information matrix is calculated as 

det Gg = {{£g,i\£gg) - 70,i)(( 4 >,2^0,2) - 70,2) 

— ((^ 0 , 11 ^ 0 , 2 ) — 70 , 170 , 2) 2 

= (^0,l|^0,l)(^0,2|^0,2) — ((£0,11^0,2) ) 2 

170,2-^0,1 — 70,1^0,2| 2 
= (m 0 |m0) - (sg x m g \sg x m e ) 

= (1 - Sg)(mg\Qgmg) 

= , J 1 2 {ng\Q g 1 ng), 

1 s g 

and hence we obtain (£ g '\Qg 1 £g) = (1 — s g ) det Gg. 


Next, we use the relation for the RLD Fisher infor¬ 
mation; g e<ij = (1 - sl)- 1 [(t e , i \t e , j ) - (2 - sD'yg^gj - 
i(£g t i\Fg£gj)] to calculate the determinant det Gg as 

(1 - s 2 g ) 2 det Gg 

= ((£0,^0,1} - (2 - a|) 7 I, 1 )«^0,2^0,2> - (2 - 4 ) 7 0 2 2 ) 

— | (£0,1 ^0,2) — (2 — 8^)70,170,2 — i(^0,i|F0£0, 2 )| 2 
= (m0|m0) - (2 - s 2 g)(sg x mg\sg x mg) - |(s0|m0)| 2 
= (! - sg) 2 {m g \Q e mg). 


This proves the relation; (1 — sg) det Gg = det Gg. □ 

2. 2y/detW \{£ g \Fg£ g ) | = TrAbs{WTmG7 *} = 

TrAbs{WTm Zg}: 

From the definition for the matrix Zg , the imaginary part 
is expressed as 


Im Zg = (£g\Fg£g) 


0 1 
-1 0 


and the straightforward calculation yields 

TrAbs {WIm. Zg} = 2Vclet W I {l\\ Fg£ 2 e ) I . 


(A4) 


(A5) 


The imaginary part of the RLD Fisher information ma¬ 
trix is 


ImG0 = 


(^0,i|-p0^0, 2 ) / 0 — 1 


1 - s 2 

1 s g 


1 0 


and the imaginary part of the inverse is 


Im{G e - 1 } = 


(^0,i|-Pg^0, 2 ) / 0 1 

det Gg l -1 0 


(A6) 


(A7) 


where we use (1 — s g ) det Gg = det Gg. It is easy to show 
(£ g \Fg£ g ) = (£gg\Fg£g t2 ) / detG 0 and thus we obtain the 
important relationship; Imx^ = ImGj 1 . This proves the 
claim. □ 

3 . mw-'jg) = det W-'Gg (Cg[W] - C?[W]) /(I - s 2 ): 

This can be shown by the following calculations: 
det W x ( 70 |VF _ 1 70 ) 

= W2270.1 + ^1170,2 - 2^1270,170,2 
— (1 ~ Sg) 1 (t02250,ll + w ligg, 22 ~ 2u>125l 2 ) 

- Re ( 102250,11 + 101150,22 - 21012512 ) 

= (1 -s^) _1 det GgTr {WGg 1 } - det GgTr jlFReG^ 1 j 
= d^ Tr { ;v . (G _, 

= ( C Z g[W ] - C?[W\) . □ 

1 s g 


3. Lemma III.3 

Since the imaginary parts of the inverse of the RLD 
Fisher information matrix and the matrix Zg are always 
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identical for two-parameter qubit models, equivalence be¬ 
tween 1 and 2 is the direct consequence of Lemma B.l. 
We next note the following relation: 


Gg 1 - ReG^ 1 


1 

det Gg 


7e, 2 -7<M70,2 \ 

-70,170,2 7e,i ) 


which follows from Eq. (61) and Lemma III.2-1. Thus 
Gg 1 — ReG ^ 1 = 0 if and only if 7 g^ = 7@ j2 = 0 . 

To show the statement about global D-invariance, we 
only need to show 70,1 = 70,2 = 0 for all 9 £ 0 if and only 
if s$ = |sg| is independent of 9. When 7 g }i = (sg\ig^) = 0 
<£> (sg\diSg) = (\/ 2 )di\sg \ 2 = 0 holds, the integration of 
this condition gives |sg| is independent of 9 = ( 9 1 ,9 2 ). 
Conversely, |sg| does not depend on 9 , we obtain 707 = 0 
(i = 1 , 2 ) for all (9 £ O. □ 


4. Lemma III.4 


When 6 = 0, the function to be minimized is /(£) = 
(£|A£) + 2|c|. Since A is positive-definite, the minimum 
is 2|c| and is attained by £* = 0. 

For the other case (6 ^ 0), introducing the new vari¬ 
ables through 


rj = BA XI % with B 


b T A "V 2 


b T A~ x ! 2 


0 -1 
1 0 


we can express the function /(£) as 

f(v) = ~(Vi +V 2 ) + 2 \Vi +c|, (A8) 

a 

where a = (b\A~ l b) is a positive constant. The minimum 
of this simple quadratic function is obtained by analyzing 
the case 7/1 + c > 0 and 71 + c < 0 separately. The result 
is 

! — a + 2c ifc > a 

— if | c| < a , (A9) 

a 1 1 

—a — 2c ifc < —a 

and the unique minimizer ?jl = ((*, 0 ) T is 

{ — a ifc > a 

—c if|c| < a ■ (A10) 

a ifc < —a 


This solution can be translated into the original variables 
£ to prove the lemma. □ 


Appendix B: Supplemental materials 


1. Canonical projection and invariant subspace 


space on real numbers. These results can be generalized 
to more general settings to be applied to the quantum 
estimation problem. 

Consider an n-dimensional vector space Lont and an 
inner product on it (•, •) : V 2 —> M. Let w±, W 2 , ■ ■ ■, Wk £ 
V (k < n) be a set of linearly independent vectors and 
define the subspace W = span R {uii, W2, ■ ■ ■, Wk} C V 
spanned by these vectors. We define a real positive semi- 
definite matrix G = an d its inverse 

by G -1 = A set of vectors jiC = 9^ w j}i= 1 

forms the dual basis of {u>i}£L 1 . 

Given a vector v £ V, the canonical projection of v 
onto the subspace W is a map kw : V —> W such that 

k 

7 tw{v) = ^^(v,w l )wi £ W. (Bl) 

i= 1 

This canonical projection is unique and it preserves the 
inner product as 

(7 tw(v),w) = (v,w) \/w £ W. Vv £ V. (B2) 

Consider any element v £ V and the condition v £ W is 
equivalently expressed as follows. 

v £ W <£> ttw(v) = v 

o Mv' £ W ± (■ v ', v) = 0 

Vt/ £ V (v 1 — v) = 0 

<=> V(u 1 ,u 2 ,..., v k ) £ V k , 

Vi(7Tw(u*) = w l => (v l — w l ,v) = 0), 

where W ± = {u £ V \ {v, w) = OVru £ W} is the (unique) 
orthogonal complement of W. 

Next, consider a linear map A from V to itself. The 
subspace is said invariant under the map A, if the image 
of W is a subset of W, i.e., A(w) £ W holds for Vru £ W. 
Using the above equivalence, this can be written as 

W is an invariant subspace under A (B3) 

<£> Vi A(wi) £ W 

ViVu' £ W 1 - (v',A(wi)) = 0 
^V(v 1 ,v 2 ,...,v k ) £ V k , 

Vi, j(7rw(^) = w l => (v l — w z ,A(wj)) = 0). 


2. Characterization of D-invariant model 

Lemma B.l Let L l g (Lg) be the SLD (RLD) dual oper¬ 
ator and Gg (Gg) be the SLD (RLD) Fisher information 
matrix, respectively. Define a k x k hermite matrix by 
Zg = [(Lg, L 3 e )+ ]i y je{i,...,k} > then, the following condi¬ 
tions are equivalent. 


In this appendix, we summarize the concept of canoni¬ 
cal projection and invariant subspace for an inner product 


1. A4 is D-invariant at 9. 

2. Vg(Ug) = ^(Im Zg^Lgj, V* £ { 1 , 2 ,..., k}. 
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3. Z e = Gg 1 

4 - £$ = £$, Vie{i,2,...,fc}. 

5. MX' 1 £ Chin), X 1 - Lg±Tg(M) with respect to 
(-, -) pe => X* — LgA-Tg(M) with respect to (-, ■)+. 

In the above lemma, {X,Y) pe := Retr (p^YX*) and 
(X, U) + := tr(peYX*) denote the SLD and RLD inner 
products, respectively. 

Proof: We prove this lemma by the chain; 1=>2=^3=^ 
4 => 5 => 1. Suppose that a given model is D-invariant, 
this is equivalent to say that the action of the commu¬ 
tation operator on the SLD dual operators is expressed 
as 

V 6 {L i e ) = Y J c ii Le , j , (B4) 

3 

with some real coefficients c*- 7 '. These coefficients are ex¬ 
pressed as 

^ = (Li,Ve{Li)) Pe , (B5) 

which directly from the orthogonality condition (13). Us¬ 
ing the relation (36), the right hand side is also expressed 
as Im Zg , and if the model is D-invariant at 0, 

Vg{U e ) = Y J ^Z e rLe, j . (B6) 

3 

Hence we show 1 => 2. 

Next, if the condition (B6) holds, the SLD inner prod¬ 
uct between Vg{Lg) and the RLD dual operator L J g gives 

(Li,Ve(L j e )) pe =J2c kj (Li, L o) P e=c ii . (B7) 

k 

The left hand side is also calculated from Eq. (37) as 
— \{g l g — g l g). Therefore, we show 

c ji = Im 4 j = -i- gl j ) ^ Gg 1 = Zg, (B8) 

holds, if the condition (B6) holds, that is, 2 =>■ 3. 

Consider an arbitrary linear operator X and assume 
the condition (B8). In this case, the RLD inner product 
between X and L l g is calculated as 

(X,t'<>)t=T'$(X,Lg, k ) + 

k 

= 5>«<X,L M )+ 

k 

= Y,4 j (X’Lg, k ) pe 

k 

= Y^(9$ j ( X ’Le,k)p»+3mz k *(X,L e , k ) P9 ) 

k 

= (X,Li) pe +i(X,Vg(Li)) pe 
= (X,Li)+, 


where Eq. (B8), and the several equations presented in 
Sec. IIC 2 are used. Since X £ C{H) is arbitrary, it 
implies 

L z e = Lg folds for all i £ { 1 , 2 ,..., k}. (B9) 

Therefore, we show 3 => 4. 

Next, let us assume the condition (B9) and we show 
that this implies the condition 5, that is, MX 1 £ 
Ch[n), X 1 — LgLTgiM) with respect to (-, ■) Pe => X* — 
LgCTgiM) with respect to (-, ■)+ . This is because 

X i -L i gYTgiM)w.r.t.{-, ■) Pe 
^Mj{X'-L l g,LgJ pe 
^Mj(X i - Li,Lg tj )+ 

^Vj(X*-4,^)+ 

^Vj(X i -4,^)+ 

«r-LJiT fi (M)w.r.t.( v )i. 

The remaining to be shown is that the condition, 

MX i £ C h {n), Mi(X 1 - UgLTgiM) w.r.t. {-,-) pe 

=► X* -UgLTgiM) w.r.t. (■,■)+), (BIO) 

implies the D-invariance of the model. Consider a 
set of hermite operators X = (X 1 , X 2 ,..., X k ) and 
suppose the above condition (BIO). Since X 1 — 
LgYTgiM) w.r.t. (-,-) Pe is equivalent to nT g (X l ) = Lg 
with 7 tt 6 the canonical projection on Tg(Xi), we can 
rewrite it as 

MX* £C h in),Mi,ji(X'-Ll,L ej ) pe =► (X i -Li,L 6iJ )+ e ). 

(BH) 

The use of Eq. (35) leads {X 1 - L^, Lg j) pe => (X* - 
Li,L ej )+ =► {X i -L} ) ,Lg ij ) P9 => iX l -n 0 ,V e iL e ^)) pe . 
Then, the equivalent condition (B3) can be applied to 
conclude that the subspace Tg)M.) is invariant under the 
action of linear operator T>g , that is, the model is D- 
invariant. □ 


3. Proof for Theorem II. 5 

i) Proof for the RLD CR bound case: 

The sufficiency (D-invariant model => C^\W] = Cg-\W\ 
for all W > 0.) follows from Proposition II.4. If the 
Holevo bound is identical to the RLD CR bound for 
all weight matrices, then all the RLD dual operators 
must be hermite, i.e., (L l 6 )* = L l g for all i = 1,2,... ,k. 
To see this, we note that X = (X 1 , X 2 ,..., X k ) £ Xg 
implies X £ Xg := {X | MiX i £ C{U), tr (pgX*) = 
0, Mi, j tr ipipgXi) = £{}. By rewriting tr (dipgXi) = 
5\ <t=> (X* — L' g ,L e ,j)+ = 0, we see that X 1 — Lg 
ii = 1,2,..., k) are orthogonal to the complex span of the 
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RLD operators, Tg(M) := spancjLe^}. It is straightfor¬ 
ward to show that Zg [X] = Zg [X — Lg] + Zg [Lg] holds for 
X G Xg, where Lg = (Lg, Lg,..., Lg) is the collection of 

the RLD dual operators. Since the matrix Zg[X — Lg] is 
positive semi-definite, and hence, we obtain 

X G Xg =4 Zg[X\ > Gg 1 2 3 4 5 6 7 8 9 , (B12) 

and the equality if and only if Vi X 1 = L l g holds. Thus, 
all Lg need to be in the set Xg. 

Next, we use the relation between the SLD and RLD 
operators, (I + i T>g)(Lg p ) = Lg^ (see Eq. (35)), to get 

(I + Wg)(Li) = Y,9 ji Le,r (B13) 

j 

Then, the conditions (Lq)* = L l e for all i imply 

Lg = y] Re{g ji }Lgj, (B14) 

j 

Vg(Li) = Y / ^{9 ji }Lg, j . (B15) 

j 

Since Gg is positive definite, so as ReGg. With these 
relations, the action of the commutation operator on the 
SLD operators are calculated as 

Vg(Lg :l ) = y^[(Re {G -1 }) -1 ]jjIm {g k ^}Lg^- (B16) 

j,k 

This proves that Vg(Lg p ) G Tg(M) = span R {Lfl ;i }. 
Therefore, if Cg\W] = C^\W] for all IT GW, then the 


model is D-invariant. 

ii) Proof for the D-invariant bound: 

This equivalence is a direct consequence of the proposi¬ 
tion II.4 (Zg = Gg 1 ) and the property of the canonical 
projection given in Appendix B 1. First, let us note 

VIT G W Cf [IT] = Gf [IT] <t=> MX G Xg Zg[X] > Zg. 

(B17) 

This is because Lg is an element of the set Xg and 
Zg[Lg] = Zg. If the condition B.l-5 holds, it is easy 
to show X G Xg => \/i,j(X l — Lg,Lgj) Pe = 0 => 
Mi,j(X 1 — L' e ,Lg t j)+ = 0. Therefore, we obtain the ma¬ 
trix inequality: 

Zg[X] = [(X\X%kj 

= [(X* - L\,Xi - LX t \ itj + [{Li,Li)+]ij 
>[(Li,Li)+] itj =Zg, (B18) 

holds for all X £ Xg because of the semi-definite positiv¬ 
ity of the matrix {(X 1 - L^X^ - L J e )+ e ] id . 

Conversely, let us assume MX G Xg Zg[X] > Zg 
is true. This is possible if and only if VX G Xg => 
Vi, j (X 1 — L l g,Lgj)+ = 0, otherwise we can always find 
some X G Xg such that the matrix inequality Zg [X] > Zg 
is violated. Thus, Lemma B.l-5 and the equivalence 
(B17) prove this theorem. □ 
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